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ABSTRACT 

Distant Type la and II supernovae (SNe) can serve as valuable probes of the history 
of the cosmic expansion and star formation, and provide important information on 
their progenitor models. At present, however, there are few observational constraints 
on the abundance of SNe at high redshifts. A major science driver for the Next Gen- 
eration Space Telescope (NGST) is the study of such very distant supernovae. In 
this paper we discuss strategies for finding and counting distant SNe by using repeat 
imaging of super-critical intermediate redshift clusters whose mass distributions are 
well-constrained via modelling of strongly-lensed features. For a variety of different 
models for the star formation history and supernova progenitors, we estimate the like- 
lihood of detecting lensed SNe as a function of their redshift. In the case of a survey 
conducted with HST, we predict a high probability of seeing a supernova in a single 
return visit with either WFPC-2 or ACS, and a much higher probability of detecting 
examples with z > 1 in the lensed case. Most events would represent magnified SNe II 
at z ~ 1, and a fraction will be more distant examples. We discuss various ways to 
classify such events using ground-based infrared photometry. We demonstrate an ap- 
plication of the method using the HST archival data and discuss the case of a possible 
event found in the rich cluster AC 114 (z = 0.31). 
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1 INTRODUCTION 

Recent observational efforts to detect fiigh-redshift (z) su- 
pernovae (SNe) have demonstrated their value as cosmologi- 
cal probes. The discovery and systematic study of faint, dis 
tant type la supernovae (SNe la) has led to rene wed progress 



in constraining the cosmic ex pansion history (Riess et al 
199S ; Perlmutter et al. 1999| ). Two independent corrected 



Hubble diagrams based on SNe la discovered via the Super- 
novae Cosmology Project (Perlmutter et al. 1999) and the 
High-Redshift Supernovae Search Team (Riess et al. 1998) 
yield a trend that strongly excludes the hitherto popular 



Such studies have now identified over 100 z > 0.2 
SNe la, and, using controlled subsets of these SNe la sam- 
ples, the first constraints are now emerging on the rate of 
their occurrence (Pain et al. 1996, 2000a). The observational 
determination of these rates is also important in many cos- 
mological applications, for example as a diagnostic of the 
cosmic star formation history (SFH) and metal enrichment 
at high- z. SNe are independent of some of the biases associ- 
ated with traditional tests based on flux-limited galaxy sam- 
ples. Such biases are increasingly of concern given the steep 
luminosity functions no w being witnessed for star-forming 



Einsteiii dc Sitter universe and, for a spatially ffaL iiiflalion- 
ary uni verse co nsistent nitli luiuiil m i crowave backgroun d- 



galax ies at all redshifts ( Steidel et al. 1999 ; Sullivan et al 



meas urements ( Hancock et al. 1998 



le Bernardis et al 



200C ), su ggests a significant non - zero cosmologica l constant, 
A ~ 0.7 (Perlmutter et al. 1999|; |Riess et al. lggj). 
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200C) and, though there remain uncertainties in the correct 
treatment of dust in high-z type II supernovae (SNe II), they 
offer an alternative approach to studying the cosmic SFH. 

SNe II are caused by the catastrophic core collapse of 
massive (> 8Mq) stars and, unlike SNe la, represent a di- 
rect probe of the massive SFH due to their short progenitor 
lifetimes - the evolution of the SNe II rate should closely 
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match tha t of the cosmic SFH ( Madau, Delia Valle & Pana- 
gia 199Kt Sadat et al. 1998). Unfortunately, due to their 



intrinsically faint, heterogeneous nature, and their likely oc- 
currence in dusty starburst regions, to date they have yet 
to be exploited in cosmological studies. However, as the pre- 
dicted SNell counts to faint survey limits (e.g. I = 26) will 
outnumber those of SNela (Madau et al. 1998; Dahlen & 
Fransson 1999, Porciani & Madau 2000), the detection of 
SNell at higher- z should become more practical as these 
fainter limit s are attained - either thr ough facilities such as 



the NGST ( |Dahlen fe Fransson 1998j ) or via studies of the 
type introduced here. 

SNela, by contrast, arise from lower mass (~ 3 — 8Mq) 
stars, and their rate depends strongly on the (as yet) un- 
known configuration of the progenitor system as well as the 
cosmic SFH. The SNela rate therefore cannot yet be used to 
directly probe the SFH without a more detailed knowledge 
of the physics of SNela explosions. The observed rate does 
however offer i mportant implications for the possibl e pro 
genitor models ( Madau, Delia Valle fe Panagia 199£ ; Yun 



gelson i jg Livio 1998[ ) and, at z > 1.5, the predicted metal 
dependency of the SNe la rate wi ll provide valuable clue s to 
the metal enrichment at high-z (Kobayahsi et al. 199£), as 
well as for the practicality of exten ding the Hubble diagram 



for more precise cosmological tests ( Ruiz-Lapuente & Canal 
199SJTI 

The study of distant SNe of both t ypes is one of the 



Frans- 



key programs underpinning the NGST (Dahlen 
son 199^). The optimum strategy for finding and studying 
SNe, particularly with respect to the investigation of SNe II 
and the SFH at z > 1, clearly depends on their abundance. 
In this paper, we suggest a strategy for producing the first 
observational constraints using existing facilities. The strat- 
egy proposed follows the idea proposed by Kolatt & Bartel- 
mann (1998) an d a similar approach used to locate faint 
sub- mm sources ( Smail, Ivison fe Blain 1997 ) . By exploiting 
the magnification bias afforded by counting sources viewed 
through the cores of strong len sing clusters, fainte r detec- 
tion thresholds can be reached (Blain et al. 1999c). In the 
case of the sub-mm surveys the intrinsic flux limits achieved 
are fainter than those reached in more ambitious dedicated 
deep surveys. Adopting the same principle, optical searches 
for lensed SNe should provide a first glimpse of the evolution 
of SNe at high-z ahead of the launch of NGST. 

A plan of the paper follows. In §2, we review the theoret- 
ical aspects of gravitational lensing as applicable to our pro- 
posed search and discuss the uncertainties associated with 
modelling SNe rates. This allows us to predict the source 
counts and redshift distribution for a survey with and with- 
out the effect of lensing. In §3 we discuss a possible appli- 
cation of the technique using the HST data archive. Only 
two rich clusters with the necessary lensing credentials have 
been multiply-visited in a single passband by HST, and 
we discuss one possible SNe found in each. In §4 we dis- 
cuss the possibilities of classifying such detections via var- 
ious follow-up strategies. We present our conclusions in §5. 
Throughout the paper, for consistency with the SNela and 
microwave background results, we assume a lambda domi- 
nated cosmology of Om = 0.3, £7a = 0.7, and h — 0.5, where 
H = 100/ikms _1 Mpc- 3 . 



2 PREDICTED SUPERNOVA COUNTS 
2.1 Introduction 

The approach proposed is identical to that succ essfully em- 



ployed in probing the faint su b-mm sources (Smail, Ivi- 



son & Blain 1997; Blain 1997). The core of a rich lens- 



ing cluster magnifies suitably-located background sources 
by factors of several. The critical difference with the sub- 
mm studies is that, for faint SNe studied at optical wave- 
lengths, the exquisite angular resolution of HST is essen- 
tial. The modest boost in detectability afforded by survey- 
ing through a massive foreground cluster significantly ex- 
tends the mean detectable redshift. We primarily consider a 
search using the HST and the Wide Field Planetary Camera 
2 (WFPC-2). The high spatial resolution, dark sky back- 
ground and relatively large field-of-view (5.71 arcmin 2 ) of 
this telescope/instrument combination ensures faint survey 
limits can be reached in a reasonable exposure time. We also 
briefly consider the advantages of using the forthcoming Ad- 
vanced Camera for Surveys (ACS) HST instrument, with a 
larger (11.33 arcmin 2 ) field-of-view. 

In the rest of this section we investigate this observa- 
tional approach by calculating the number density N(m, z) 
of SNe of different types with and without the magnification 
produced by a massive cluster located at various redshifts. 
Below we introduce the theoretical aspects of both the lens- 
ing and the intrinsic SNe properties as a function of redshift. 



2.2 Cluster Mass Models 

The role of clusters in gravitational lensing has been re- 
viewed b y ma ny authors (see, for example, Blan dford & 
Narayan ( 1992- ) or Hattori, Kneib & Makino ( |l999| ) and ref- 
erences therein). A gravitational lens will distort and mag- 
nify a background source by an amount that depends on 
both the position of the source behind the cluster, and the 
relative redshifts of the source and the lens. This results 
in different background regions suffering different degrees of 
amplification, depending on their location behind the clus- 
ter. Through the locatio n and spectro scopic identifications 
of multiply-imaged arcs (Mellier 1999), it is possible to cor- 
rect for this effect by utilising a well-constrained cluster mass 
model to the accuracy required. These models involve multi- 
ple component mass distributions that are able to trace the 
potential well of a cluster and its m ore massive individual 
member galaxies (Kneib et al. 1996). The unce rtainties in 



these best cases is currently ~ 10 — 20 per cent (Blain et al 



1999a) 



The degree of magnification that one can expect in the 
central regions of such massive clusters - the areas sampled 
when conducting lensed SNe se arches - are typical ly of order 
H ~ 3-4 for sources at z s — 2 (Blain et al. 1999a), where fj, 
is defined as the ratio of the magnified flux to the flux mea- 
sured in the absence of lensing. However, this boost in sen- 
sitivity - or effective increase in the SNe number density as 
fainter SNe become brighter than the survey limiting magni- 
tude due to amplification - is offset by the competing effect 
of depletion. This results in a decrease in the SNe number 
density, as the unit solid angle in the unlensed source plane 
increases by a factor of l//i(z) in the lensed case, resulting 
in a decrease in the SNe number density by l//i(z). 
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When attempting to determine a precise count rate 
of lensed sources, accurate modelling is essential. Here, we 
choose to model the cluster mass distribution more sim- 
ply via a singular isothermal sphere (SIS), where the mass 
density p of the virialized dark matter halo of a cluster is 
given by p(r) — a 2 /2ivGr 2 , where r is a radial coordinate, 
and p is parameterized through a v , the one-dimensional 
velocity dispersion. The mass M of a halo interior to r 
is then M(r) — 2a 2 r/G. This corresponds to a constant 
deflection angle for incident light rays, /3 = 4ir(a v /c) 2 — 
28.8" (a\,/1000 km s -1 ) 2 , always pointing towards the lens 
centre of symmetry. For a source that, in the absence of lens- 
ing, would be seen at an angular distance 8 from this centre, 
the lens equation leads to an image at 8+ = 8 + 8e, with 
magnification p+ = Qe/8 + 1, where 8e = PDi s /D oa is the 
Einstein radius, and D oa and Di 3 are the angular diameter 
distances between the observer, the lens, and the source. If 
the alignment is 8 8e (strong lensing), a second image is 
produced at 8- = 8 — 8e with magnification /i_ = Oe/0 — 1. 

High-resolution iV-body simulations have recently 
shown, however, that haloes formed through hierarchical 
clustering have a universal density profile which is shallower 
than isothermal (but still diverging like p oc near the 

halo centre, and steeper than isothermal (with p oc r~ 3 ) in 
its outer regions (Navarro, Frenk, & White 1997, hereafter 
NFW). This profile takes the form: 

p c 5 C 



p{r) = 



(r/r s )(l +r/r a ) 2 



(1) 



where p c = 3H (z)/87rG is the critical density at the lens 
redshift, <5 C is a characteristic over-density, and r s = r2oo/c is 
a scale radius, where c is a dimensionless number measuring 
the "concentration" of the halo, and the virial radius, r2oo, is 
defined as the radius inside which the mass density is equal 
to 200p c . S c is related to c via 



200 



3 ln(l + c) 



(2) 



l+c 



We can then see that M(r 20 o) = 8007rp c r-| 00 /3. The NFW 
haloes are thus defined by two parameters; c, and either r2oo 
or M 2 oo- 

The lens equation for the NFW density distribution has 



been c omputed by several authors (Bartelmann 1996 ; Maoz 
et al. 1997; Wright & Brainerd 1999), and a go od su mmary 
of its application can be found in Hattori et al. (1999). With 
respect to SIS profiles containing the same total mass and 
located at the same redshift, NFW lenses have smaller cross 
sections for producing multiple images, but this is compen- 
sated for by the fact that the average magnification of the 
images is higher for NFW haloes. In other words, NFW 
haloes are efficient magnifiers but poor image-splitters - a 
SIS profile underestimates the magnification obtained com- 
pared to an NFW profile. 

In the following calculations, we model the lensing clus- 
ters using both SIS and NFW profiles. As a guide to the 
benefits of searching for SNe through clusters, we have ex- 
plored two hypothetical cases: a survey conducted through 
a sample of (known) a v = 1000 kms -1 z ~ 0.2 clusters, 
and a similar sample of z ~ 0.6 clusters now being located 
through the MACS Survey (see for example Ebeling, Edge 
& Henry 2000). To meaningfully compare the magnifications 
of SIS and NFW profiles, we assume that for both profiles 




Point source (SNe) redshift 

Figure 1. The median area- weighted magnification suffered by 
a point source as a function of redshift behind typical a v = 
1000 kms -1 , z = 0.2 and z = 0.6 clusters, for a search conducted 
using WFPC-2. The amplifications expected for both NFW and 
SIS profiles are shown. 



the clusters have identical masses interior to the virial ra- 
dius, T2oo- We then equate the expressions for M(r) at r2oo 
for both profiles, and hence set r2oo for the NFW profile 
given a value for a v . For the NFW profiles, the concentra- 
tion parameter c for each cluster is calculated for our cos- 
mology using the FORTRAN routine charden.f provided by 
J. Navarro. 

The median area- weighted amplifications are calculated 
by considering a pixel matrix corresponding to the WFPC-2 
CCD, and determining the amplification suffered by a back- 
ground source were it to appear in each element. The median 
pixel amplification can then be found. Fig. jj] shows this am- 
plification of a background source as a function of source 
redshift for two typical clusters. Typical median amplifica- 
tions of p ~ 2-3 are typical of z ^ 1.5 sources, corresponding 
to a brightening in magnitudes of ~ 1. 



2.3 Supernova Rates 

We now estimate the number of SNe we would expect to find 
to a given survey magnitude limit at a given epoch, in the 
case of a blank field (i.e. without any significant lensing mag- 
nification) and through a rich cluster. We will not include 
those (unlensed) SNe that occur in the cluster galaxies, and 
will therefore underestimate the total number of new SNe 
that will be found in a lensed search. The number counts 
depend on both the intrinsic redshift-dependent SNe rate, 
and that period of the light curve during which time the 
SNe appears brighter than the survey search limit. 

A number of workers have modelled the expected evolu- 
tion in the SNe rate with redshift. Most have followed a sim- 



ilar route to that employed by Madau et al. (1998), who base 



their estimates on the observed luminosity densities which 
trace the cosmic star formation history (SFH). Using various 
diagnostics, the cosmic SFH can now be traced to high red- 
shift [z — 4), although some details remain controversial. As 
the progenitor lifetimes for SNe II are very short by cosmo- 
logical standards, the SNe II rate can be calculated directly 
from the SFH provided an initial mass function (IMF) is 
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assumed. Here we adopt a Salpeter (1955) IMF with lower 



and upper mass cut-offs of 0.1 Mq and 125 Mq respectively. 

In modelling the SNe II rate, the form of the cosmic SFH 
is critical. The original study of the evolution in the comov- 
ing star- formation rate base d on optically selected galaxy 
surveys ( [Madau et al. 1996 ) suggested a sharp rise from 
z = to a peak at z ~ 1-2 by a factor of almost 1 (Lilly et 
al. 199f]; [Madau et al. 1996| ; |Connolly et al. 1997| ), followed 
by a decline at high-z. However, more recent studies suggest 
that the evolution of the SFH up to z ~ 1.5 may have been 
overestimated, perhaps only increasing by a factor of ~ 4 



(iTrcssc fc Maddox 199^ : 



Cowie, Songaila & Barger 1999 



Sullivan et al. 200C), and that the SFR at high-z may be 



obscured due to large amounts of dust extinction (see, for 
example, Pettini et al. 1998 and Blain et al. 1999b). The 
form of this evolution to z ~ 2 is crucial for our studies 
when predicting the SNe II rate. 

We consider two illustrati ve SF Hs. The first (SFH-I) is 
taken from Madau & Pozetti (2000), and provides a conve- 
nient analytical fit of the form: 



Pb{z) = 



0.23 e J 



e 3.fe +44.7' 



(3) 



which assumes an IMF lower mass cut-off of ~ 0.5 Mq and 
an Einstein de Sitter (EdS) Universe. Multiplying by a fac- 
tor of 1.67 will convert to a Salpeter IMF with a lower mass 
cut-off of 0.1 Mq. We convert the SFH to the A-dominated 
universe assumed here by computing the difference in lumi- 
nosity density between an EdS Universe and our cosmology, 
and applying this correction to the fit above. The SFH fit 
matches most UV-continuum and Ha luminosity densities 
from z — to z — 4, includes a mild correction for dust of 
^4i500 = 1.2 mag (A2800 = 0.55 mag), and implies a local 
SNe II rate of 6.2 x 10~ 5 SNeyr" 1 Mpc 3 , v ery cl ose to that 
observed by Cappellaro, Evans & Turatto ( 1999 ). 

However, the SFR evolution in this model to z ~ 1.5 
is both stronger, and results in a lower local SFR, than 
more recent measurements (Cowie et al. 1999; Sullivan et 
al. 2000). Accordingly, we also consider a second SFH (SFH- 
II) with a shallower evolution (a factor of ~ 4 from z — to 
z ~ 1.75 in an EdS Universe, and constant thereafter), but 
one which over-produces the local SNe II rate. The SNe II 
rates as a function of redshift are shown in Fig. ^. 

The SNe la rate is more difficult to model due to the 
uncertain nature of the progenitor system. The evolution 
expected depends critically on whether SNe la occur in dou- 
ble or single de genera te progenitor systems (for a review see 
Nomoto et al. (1999) and references therein), and there is 
much debate in the literature as to what form any evolution 
should take. We consider a commonly used empirical ap- 
proach, and represent the uncertainty in terms of two param- 
eters (Madau et al. 1998; Dahlen & Fransson 1999). The first 
is a delay time t, between the binary system formation and 
SNe explosion epochs, which defines a (time-independent) 
explosion probability per white dwarf. The second is an ex- 



plosion efficiency, 77, which accounts for the fraction of bi- 



nary systems that never result in a SNe. This efficiency r\ is 
constrained using a preliminary determination of the SNe la 
rate a t z = 0.55 of 1.66 ± 0.42 fe 3 10~ 4 Mpc~ 3 yr~ 1 by Pain 



SNe II rates 



/ 

/ 



r=0.3 Gyr 



O L<*' /'/ 



T=3.0 Gyr 



1 2 

Redshift 



Figure 2. The adopted SNe rates as a function of redshift for 
both SNe II (top) and SNe la (bottom). The long dashed lines 
refer to rates derived using SFH-I, the short dashed lines refer 
to SFH-II. The SNe la rates are shown using two 'delay' times 
(t = 0.3 & 3.0 Gyr). The horizontal solid line and shaded zone 
indicates the z = 0.55 SNe la rate and its uncertainty obtained by 
Pain et al. (2000a), used to calculate the efficiencies of the SNe la 
rates. Local SNe rates determined by Cappellaro et al. (1999) are 
shown for comparison. 



first has a small value of t = 0.3 Gyr (corresponding to a 
shallower decline at high-z,) whilst the second is a larger 
value of r = 3.0 Gyr, which produces a steepe r drop -off. For 
comparison, the z = 0.55 rate of Pain et al. (2000a) is also 
shown in Fig. |^, together with the difference that arises from 
using the two SFHs, and how these histories tend to over or 
under-predict the local SNe rates. 

While much published data exists for SNe la, by con- 
trast there is a dearth of data for SNe II and a great variety 
in their spectra, peak magnitudes and light curves. For the 
purposes of this p aper , we adopt the models of Gilliland, Nu- 
gent & Phillips (19991) (hereafter GNP) which incorporate 
recent light curve templates, fc-corrections and luminosity 
functions for SNe la and II. These models include observed 
dispersions in the peak magnitudes of SNe la and II, which 
are included in our simulations of our SNe counts (see GNP 
for a full discussion of the models and simulations). 

The major advantage to using GNP's models lies in the 
fact that the Spectral Energy Distributions (SEDs) are more 
realistic than those employed previously for rate calculations 
and follow the true temporal evolution of these objects. The 
SNe la template was constructed using data from IUE and 



HST (Cappellaro, Turatto & Fernley 1995 



Kirshner et al 



1993) in the UV and a large number of ground-based ob- 
servations in the optical. Spectrum synthesis calculations 
were carried out by P N for the SNeH's using t he non-LTE 
code PHOENIX 9.0 ( |Hauschildt et al . 1997|; Irlauschildt ' 



et al. (|2 | 000a|) using a sample of 38 SNe la from the SCP (see Baron, fc Allard 1997|). The s pectra of SNe 1979C ( |Cap~- 



also Pa n et al. 2000b). 



laro, Turatto fc Fernley 1995), 1992H and 1993W (courtesy 



We consider two SNe la evolutions, shown in Fig. The D. Leonard) were fit at several epochs. The resultant fits of 
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these spectra provided a nice match to the observed data 
and yield a crude SED template for the SNell's. The major 
emphasis of this effort was to model the UV as well as pos- 
sible (this portion of the SED dominates in high-z searches 
using filters bluer than lfim), since the simple assumption 
of blaekbody SEDs has been shown to be quite inaccurate 
(see fig. 10 of GNP). 



2.4 Results 



We can now combine the two m ass m odels from §2.2 and the 



various SNe rate models from £2.3 to estimate the number 



counts of SNe with and without the presence of a strong- 
lensing cluster. 

As with blank-field surveys, searches such as those out- 
lined here require two images of the same rich cluster, but 
with the two epochs of observation separated by a suitable 
amount to allow the SNe to become visible. When estimat- 
ing the simple count rate of lensed high-z SNe, therefore, 
both forward (SNe in second image) and backward (SNe in 
first image) searches can be included, and to allow for the 
time dilation of the high-z SNe, we assume a time of at least 
2yrs between the two exposures. Table ^ gives the counts 
expected to survey limits of 7gi4 = 26.0 & 27.0, assuming 
Zciustcr = 0.2 & 0.6. 

When comparing lensed and unlensed counts, one of 
the key factors is the slope of the SNe number counts, 
q = dlogN/dm. For a > 0.4, we expect lensing to in- 
crease the number counts, while for a < 0.4, lensing will 
have the opposite effect (depletion). The slopes in our mod- 
els depend on i) the intrinsic SNe properties, contained in 
the models of GNP, and ii) our assumed redshift evolution 
of the SNe number density. The models of GNP include the 
best-determined uncertainties in the SNe properties, for ex- 
ample in the dispersion of the peak SNe magnitude, and our 
simulations therefore account for this first uncertainty. The 
redshift evolution of the SNe numbers depends on which 
SFH we use, and, for SNe la, the value of the delay time, r. 
By running simulations for the unlensed case, we find that 
typical values of a for SNe II for SFH-I are a ~ 0.35, and 
for SFH-II q ~ 0.32. There is little variation in these values 
between simulation runs. 

The number- redshift relation is plotted in Fig. |j| As 
one would expect, the figure shows there is little difference 
in the rate of SNe up to z c i us t cr . From z c i U stcr to z ~ 1, 
the lensed counts are lower than those appropriate for blank 
field searches. This is due to the effect of depletion men- 
tioned above, which unfortunately negates the lensing ben- 
efits. However, the important differences begin at z > 1, 
where we sample a fainter, more numerous population of 
SNe. Lensing has the principal effect of increasing the pro- 
portion of SNe found at high-z, bringing fainter, magnified 
examples into view. There is then a corresponding extension 
in the redshift distribution. For example, for the z = 0.2 
NFW cluster, the mean redshift of detection changes from 
z ~ 0.70 to z ~ 0.85, and the number of detections at z > 1 
increases from 20 to 31 per cent. 

The numbers drop sharply by z = 2, even in the lensed 
case, due to the UV drop-off in the SNe spectrum being 
redshifted into the Isuw filter, which prevents the detection 
of higher-z SNe. We note that the use of HST and an infra- 
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Figure 3. Redshift distribution of a search for lensed SNe gener- 
ated from a Monte-Carlo simulation using the models of §2.2&2.3, 
to the limit J 814 = 26.0. Both SNe II (heavy lines) and SNe la 
(light lines), magnified (solid line) and unmagnified (dashed line), 
are shown. The position of the (NFW profile) lensing cluster is 
shown by the dashed line at z = 0.2. Lensing alters the redshift 
distribution, resulting in the detection of a greater proportion of 
high-redshift SNe of both types. Only forward searches are in- 
cluded here. 



red search would allow the detection of lensed z > 2 SNe 
which would be undetectable in the unlensed case. 



3 APPLICATION TO THE HST ARCHIVE 

Clearly, the number of detectable SNe depends on many 
variables which forms a strong motivation for this search, 
particularly in view of the need to quantify the high red- 
shift rate in anticipation of detailed studies with NGST. A 
robust prediction from Table [j] is the high chance of discov- 
ering a new SNe during a repeat exposure of a previously 
studied lensing cluster. Specifically, to a reasonably modest 
magnitude limit (achievable in ~ 1-2 orbits), we can expect 
to discover 1 to 3 SNe per visit if we combine forward and 
backward searches. 

To demonstrate the potential of such a search, we have 
analysed the HST data archive in the case of repeated ob- 
servations of strong-lensing clusters for which mass models 
are available. Though there are few such repeated visits with 
the same filters, exposure times, enough images for a thor- 
ough cosmic-ray rejection and a satisfactory delay between 
the 'reference' and 'discovery' images, we have none the less 
discovered a candidate event in the rich cluster AC 114. The 
core of cluster AC 114 (z = 0.31) was observed on three 
occasions with WFPC-2 in the F702W filter (see Natara- 
jan et al. 1998). The first two observations have total ex- 
posure times of 16.8 ks each and were taken on JD 50086.3 
and 50089.3 (GO#5935). Both of these observations suffered 
from scheduling problems which resulted in the images hav- 
ing ~ 4-5 x the nominal sky background, limiting their sen- 
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Table 1. The predicted number of SNe to survey limits of I&14 = 26.0 & 27.0 per WFPC-2 field for both a blank field and for a 
lensed survey (assuming zi cns = 0.2/0.6 and the two lensing models). In each column the main entries refer to SFH-I and the values in 
parenthesis to SFH-II. The totals are the sum of IIL, IIP and a mean of the two la values. The bottom two lines show the expected 
counts for a survey using ACS and one of the scenarios. Also shown are the percentage of SNe found at z > 1.0 & z > 1.5. Both forward 
and backward searches are included. 



Survey 


Lensing 






SNe la 




SNe II 




Totals 


Fraction (%) of SNe II with: 


Limit 


Model 




T - 


-- 0.3 


T = 


3.0 


(IIL+IIP) 






z > 1.0 


z > 1.5 


I = 26.0, WFPC-2 


Unlensed 





32 


(0.25) 


0.22 


(0.23) 


1.44 (1.68) 


1 


71 (1.92) 


23.2 (14.4) 


1.4 (1.2) 




SIS, z = 0.2 





25 


(0.19) 


0.18 


(0.18) 


1.30 (1.44) 


1 


52 (1.63) 


30.3 (20.1) 


4.0 (3.3) 




SIS, z = 0.6 





31 


(0.24) 


0.22 


(0.22) 


1.47 (1.68) 


1 


74 (1.91) 


26.4 (16.1) 


2.6 (1.8) 




NFW, z = 0.2 





23 


(0.18) 


0.16 


(0.17) 


1.25 (1.40) 


1 


45 (1.58) 


31.0 (20.8) 


4.6 (3.8) 




NFW, z = 0.6 





.30 


(0.23) 


0.21 


(0.21) 


1.46 (1.69) 


1 


72 (1.91) 


26.9 (17.5) 


3.2 (2.5) 


I = 27.0, WFPC-2 


Unlensed 





56 


(0.41) 


0.38 


(0.37) 


2.86 (3.10) 


3 


33 (3.49) 


32.7 (22.6) 


3.8 (2.8) 




SIS, z = 0.2 





42 


(0.30) 


0.27 


(0.27) 


2.34 (2.42) 


2 


69 (2.71) 


37.7 (26.8) 


6.4 (5.9) 




SIS, z = 0.6 





51 


(0.38) 


0.36 


(0.34) 


2.85 (3.05) 


3 


29 (3.41) 


34.5 (23.2) 


5.4 (4.2) 




NFW, z = 0.2 





.38 


(0.29) 


0.25 


(0.25) 


2.20 (2.31) 


2 


52 (2.58) 


37.9 (27.1) 


7.1 (6.3) 




NFW, z = 0.6 





49 


(0.37) 


0.33 


(0.33) 


2.77 (2.97) 


3 


18 (3.34) 


33.5 (23.4) 


5.4 (4.7) 


I = 26.0, ACS 


Unlensed 





63 


(0.50) 


0.44 


(0.46) 


2.86 (3.33) 


3 


39 (3.83) 


23.2 (14.4) 


1.4 (1.2) 




SIS, z = 0.2 




- (0.40) 




(-) 


- (2.92) 




— (3.32) 


- (19-7) 


- (3.0) 



sitivity to -R702 ~ 26 rather than the expected -R702 ~ 27. 
As a result the cluster core was re-observed in Director's 
discretionary time (GO#720T) on JD 50747.5 for a further 
15.6 ks. 

Each of the observations comprises 6 single-orbit ex- 
posures in 3 pairs, each of which was spatially offset by 3 
WFC pixels to enable the removal of hot pixels and sim- 
ilar artefacts. The observations are shifted so as to place 
the cluster cD on WFC4 (JD 50086.3), WFC2 (JD 50089.3) 
and on the boundary of WFC3-WFC4 (JD 50747.5) with 
space-craft orientations of PA(V3)=270, 270 and 220 de- 
grees respectively. The three observations of the cluster core 
therefore sample very different areas of the WFPC-2 detec- 
tors, and the combination of this with the staggered timing 
makes this dataset well suited for searching for faint variable 
sources. 

The data was reduced using the standard STScI 
pipeline, aligned with integer pixel shifts and then stacked 
using the STSDAS task CRREJ. The observations were cal- 
ibrated usi ng th e F702W zero-point quoted by Holtz- 
man et al. (1995). Finally, the three observations were ro- 
tated and aligned to allow a search for objects whose bright- 
ness varies between the visits. 

The overlapping area between the three observations 
comprises a total area of roughly 80" x 80" and this re- 
gion was searched for sources which were seen in both the 
JD 50086.3 and 50089.3 observations, but which were miss- 
ing from the JD 50747.5 observation. The sequence is shown 
in Fig. ^J. A variable point-source was identified - 4.6" 
(45kpc) out in the halo of the central galaxy. The source 
brightness is R702 = 26.4 ± 0.15 and 26.2 ± 0.15 in the ear- 
lier observations, but is fainter than -R702 > 27.6 (3a) in the 
observations taken 21 months later. 

Though the object could be a SNe la in the envelope of 
the central cD galaxy, this would make it very faint unless 
caught at the end of its light curve. Note that the absence 
of a counterimage is not unexpected given the likely time 
delay involved. 

Gal- Yam & Maoz (2000) have conducted a similar 
comparison of the strong-lensing cluster A 2218 which re- 



vealed another potential SNe. They detect a transient ob- 
ject (Veoew = 22.8) in an F606W March 1999 image, but 
which is absent from both an F606W January 2000 and an 
F702W September 1994 image. Though the relative bright- 
ness of this SNe suggests that it occurred in a z = 0.1753 
cluster member galaxy, the usefulness of strategically sur- 
veying strong lensing clusters, for both lensed and unlensed 
SNe, is well demonstrated. 



4 DISCUSSION 



As we have seen in §2.4, there is a high degree of likelihood 
of detecting a new SNe in a repeat visit to a strong lensing 
cluster. We will now discuss possible follow-up strategies to 
any detections, and the type of science that can be achieved 
from these types of surveys. 

As with earl y ground-based SNe la searches conducted 
in rich clusters ( Hansen et al. 1989| ), it could be difficult 
to isolate the SNe, constrain its type, and determine a red- 
shift. These early SNe searches were, however, essential in 
motivating astronomers to execute the later more ambitious 
programs, and before any concerted follow-up programs can 
be developed, it is first necessary to determine how many 
lensed SNe there may be and, where possible, their proper- 
ties. 

One of the primary benefits of the lensed search is the 
fact that faint, high- 2: SNe - which would otherwise be too 
faint for spectroscopic follow-up - are brightened sufficiently 
for spectroscopy to become feasible. We estimate that ~ 40 
per cent of z > 1.5 sources could be followed in this way, 
compared to ~ 5 per cent without lensing. Such highly- 
lensed SNe such as these will typically be identifiable via 
their proximity to the cluster core. However, the magnitude 
of the faintest detected objects is still too faint for spectro- 
scopic follow-ups on existing equipment. We need to exam- 
ine various other methods of exploiting discoveries of lensed 
high-z SNe, and in particular, the possibilities of distinguish- 
ing the type of SNe in order to preferentially select either 
SNe la or II. 
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Figure 4. A demonstration of the technique applied to HST archival data on the lensing cluster AC 114 (z = 0.31). The detection 
R702 = 26.3 is found in both reference images on different parts of the detector (and hence cannot be an instrumental defect), but it is 
invisible (-R702 > 27.6) in the image taken 21 months later. 



New SNe will be caught at various stages on their light 
curve. We consider a statistical approach to classifying the 
detected SNe based on their colour distribution using a 
ground-based follow-up to any Jgi4 detections using, as an 
example, the NIRC on the Keck-II telescope. Fig. H shows 
the (/ — K) - redshift distribution of detected SNe in the 
HST search using a simulation of many repeat visits to a 
2 — 0.2 cluster, with the search criteria described in §2.4 
The 'ridgeline' of blue SNe la colours arises from those nearer 
maximum light, with the SNe lis redder than the la's. The 
ridgeline of redder SNe lis arises from the fact that SNe lis 
settle down to a constant colour after ~ 50 days in the rest- 
frame. The other ridges and gaps in the diagram are a prod- 
uct of the light curve, for example there is a greater likeli- 
hood of detecting a SNe in a plateau phase than during a 
steep decline. 

Fig. |^ indicates that the bulk of the reddest objects 
would be expected to be SNe with 2 > 1. These red objects 
would be detected in follow-ups using large ground-based 
telescopes in the Jf-band; i.e. those detected in K are likely 
to be higher- 2 and likely to be SNe II. Several photometric 
measurements can determine the rate of decline over a lu- 
nation and further constrain the type. The bluest detections 
will be SNe la observed close to maximum light. These may 
be further distinguished via the presence of a host galaxy 
whose redshift can be measured. 

Although simply determining the number counts in 
searches such as these is important, it is interesting to con- 
sider other applications and their impact on our understand- 
ing of SNe. 

The number counts of SNe la predicted is small, and it 
would be difficult to advance the existing Hubble diagram 
using this method. Accurate amplification corrections to the 
detected magnitudes would also be required. The mass mod- 
els of most rich clusters are currently accurate to around 10- 
20 per cent, which would lead to uncertainties in the mag- 
nitudes of around 0.5 mag, insufficient for distinctions of a 
cosmological world model to be made. However, the mere 
detection of a high- 2 (and therefore red) SNe la would have 
important implications for the metallicity of the progenitor 
system. In a survey conducted to I — 26.0 (27.0) across 40 
clusters, we would expect 5(8) SNe la to be found in the for- 
ward search, 2-3 of which would be at z > 1. Using these, 
the first tests of the high-z SNe la rate could begin. 

A significant achievement would be the determination 
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Figure 5. The simulated I — K colour-redshift distribution for 
likely detections of SNe in a forward search behind a z = 0.2 rich 
cluster to Isi4 ^ 26.0. Each point represents a different detection, 
taking account of visibility, fc-corrections, dust (SNe II only), and 
lens magnification effects. The lensing assumes the NFW pro- 
file discussed earlier. The distribution indicates the feasibility of 
ground-based follow-ups. 



of the high- z SNe II rate, no measurement of which yet 
exists. The SNe la rate study of Pain et al. (2000a) at 
z = 0.55 shows that a sample of 38 SNe results in statis- 
tical and systematic uncertainties in the SNe rate of ~ 20 



per cent. Again, in a survey of 40 clusters, we would expect 
around 35(60) SNe II to be detected in the 'forward' search 
which could then be followed up from the ground. Of these 
~ 10(17) would be SNe II at z > 1, and using the cluster 
mass models and an assumed cosmology, measures of the 
SFR at 2 = 1-2 could be made in an independent manner 
to existing measurements. 

All these rates will be increased by using a larger field- 
of-view instrument, for example the 11.33 arcmin 2 Advanced 
Camera for Surveys (ACS). As an example, we repeated the 
simulations using this instrument and the r = 0.3 Gyr, SFH- 
II, / = 26.0 survey limit scenarios, shown as the bottom line 
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Figure 6. The distribution in I — K of detected SNe for both la 
and II types. The shaded regions indicate the colours of the high- 
est redshift objects (z > 1) and hence how follow-up strategics 
designed only to select high-z candidates might be imagined. 

in Table [l]. With this instrument, the predicted counts to 
I — 26.0 are as high as using the WFPC-2 to a survey limit 
of I — 27.0, and it becomes almost certain that a SNe will 
be found, and the larger field-of-view, and hence coverage 
of the lensing cluster, results in the potential discoverey of 
more high-z SNe. 



(v) We briefly consider the merits of using the Advanced 
Camera for Surveys, and find that even to I = 26.0, we 
expect ~ 3 SNe per search (one of which will lie at z > 1) 
with a greater extension in the redshift range probed. 



NOTE ADDED IN PROOF 

We note the preprint of Saini et al. (2000), which also ex- 
amines the effect of gravitational lensing by massive clusters 
on high-redshift SNe, in particular those located near to gi- 
ant arcs where the lensing effect is strongest. For the case of 
the cluster A 2218, they conclude such arcs could enhance 
the S/N of detected SNe by ~ 15, and that such lensing 
would substantially extend the maximum redshift of SNe 
detection. We note that they also advocate the monitoring 
of such systems for the detection of cosmologically useful 
SNe at z ~ 1 and beyond. 
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5 CONCLUSIONS 

In this paper, we have presented a technique for detecting 
SNe at z > 1 by utilising the magnification bias obtained by 
searching through rich clusters of galaxies. We summarise 
our main conclusions as follows: 

(i) We have shown that a search utilising repeat HST 
WFPC-2 visits covering the cores of rich z = 0.2-0.6 clusters 
would reveal ~ 1-3 new SNe, of which ~ 80 per cent will be 
of type II. Such a search will both increase the mean redshift 
of detected SNe, and will result in a threefold increase in the 
number of z > 1.5 SNe detected. 

(ii) We have applied our technique using HST archival 
data, and have discovered a possible SNe event in the cluster 
AC 114 at z = 0.31. 

(iii) We have discussed various follow-up strategies. The 
brightest events (to / = 24) can be followed spectroscop- 
ically using large ground-based telescopes, and a fraction 
of these will be highly lensed z > 1.5 SNe, which it would 
not be possible to follow-up without the lensing boost. For 
fainter SNe, A"-band imaging will detect and select the red- 
der, high-z SNe, and the light curve measurements can con- 
strain the type. 

(iv) We have discussed the potential scientific benefits of 
such a search, including the determination of the high-z SNe 
rates, constraints on the cosmic SFH, and the possibility of 
eliminating some SNe la progenitor models. 
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